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ABSTRACT 

We present a study of satellites in orbit around a high-resolution, smoothed parti- 
cle hydrodynamics (SPH) galaxy simulated in a cosmological context. The simulated 
galaxy is approximately the same mass as the Milky Way. The cumulative number 
of luminous satellites at z = is similar to the observed system of satellites orbiting 
the Milky Way although an analysis of the satellite mass function reveals an order of 
magnitude more dark satellites than luminous. Some of the dark subhalos are more 
massive than some of the luminous subhalos at z = 0. What separates luminous and 
dark subhalos is not their mass at z = 0, but the maximum mass the subhalos ever 
achieve. We study the effect of four mass- loss mechanisms on the subhalos: ultraviolet 
(UV) ionising radiation, ram pressure stripping, tidal stripping, and stellar feedback, 
and compare the impact of each of these four mechanisms on the satellites. In the 
lowest mass subhalos, UV is responsible for the majority of the baryonic mass loss. 
Ram pressure stripping removes whatever mass remains from the low mass satellites. 
More massive subhalos have deeper potential wells and retain more mass during reion- 
isation. However, as satellites pass near the centre of the main halo, tidal forces cause 
significant mass loss from satellites of all masses. Satellites that are tidally stripped 
from the outside can account for the luminous satellites that are lower mass than some 
of the dark satellites. Stellar feedback has the greatest impact on medium mass satel- 
lites that had formed stars, but lost all their gas by z = 0. Our results demonstrate 
that the missing satellite problem is not an intractable issue with the cold dark matter 
cosmology, but is rather a manifestation of baryonic processes. 

Key words: galaxies: dwarf — cosmology: theory — galaxies: evolution — methods: 
N-Body simulations — methods: numerical 



1 INTRODUCTION 

The A Cold Dark Matter (ACDM) cosmology is currently 
the most widely accep ted and successful paradigm for de- 
scrib i ng the Universe jBIumenthal et aT [ll984l: Id avis et al.l 
1 19851 : lGramann|[l98i : IPeebles fc Ratrall2003l '). Consequently, 
structure in the Universe forms hi erarchically, as shown an- 
alytically and in s imulations (e.g. P ress fc Schechteilll974l : 
I White fc Re3[l978 : Davi s et al.. .19851 ). First, dark matter 
collapses into small halos, and later these collect as subha- 
los into galaxies, where gas cools into a disk to form stars. 
In especially dense regions of the Universe galaxies bind to- 
gether gravitationally into galaxy clusters. 

In a hierarchical Universe, substructure is expected to 
be invariant at all scales of interest. In some of the ear- 
liest simulations that were able to resolve substructures. 



I Moore et al.l (| 19991 ) found that dark matter-only simulations 
of galaxies and galajcy clusters had the same number of sub- 
structures relative to the total mass of the system. A com- 
parison of the simulations to observations showed that the 
simulated galaxy cluster matched the quantity of substruc- 
ture in the nearby Virgo cluster, but that the simulated 
galaxy had significantly more substructure than the Local 
Group. Using con strained simulations of a system similar to 
the Local Group, iKIvpin et al.l l| 19991 ) found far more sub- 
structure than what has been observed. This discrepancy 
between ACDM and observations is kn own as the "miss- 
ing satellites problem" . .Kravtsovi (|2010l ) provides a recent 
review of the progress made towards solving this problem. 

Large observational surveys ha ve also discov e red a 
new class of ultra-faint galaxies (IWillman et al.l l2005l : 
iBelokurov et al.ll2007l : IKoposov et al.ll2008l ). The detection 
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of such galaxies slightly lessens the number of satellites that 
are "missing" from the Lo cal Group. Early obs ervations of 
their velocity dispersions (jSimon fc GehalfioOTi ) show that 
stars may form in halos of lower mass than previously be- 
lieved possible, although the typical star formation effi- 
ciency in these low mass halos must be extremely low given 
that the halo mass function rises stee ply at these masses 
l|Tollerud et all l2009l : ICuo et alJboiC f). A key question is 
whether there is a minimum mass halo in which stars can 
form, and if there is significant scatter in the star formation 
efficiency at a given halo mass. The low mean star formation 
efficiency at these masses might be driven by halo-to-halo 
variation, or a steep, universal relation between star forma- 
tion efficiency and halo mass. The ultra faint satellites have 
been detected down to My ~ —2, which corresponds to 100 
Lq , below the resolution of the simulations studied here, ft 
is impossible to give a full census of such objects from the 
simulations. However, we show that some small objects do 
form in the simulations. 

Generally, there are two paths pursued to solve 
the missing satellites problem. One is to alter the 
cosmological paradigm. Examples of this include self- 
interacting dark matter i n which subhalos are dest royed 
through self-annihilation l|Spergel fc Steinhardtl [2OOOI ) . ini- 
tially warm dark m atter out of which sma ll struc- 
tures do not form (jDalcanton fc HoganI l200ll ). or re- 
moving small s cale perturbations from the primordial 
power spectrum (IZentner fc Bullockll2003l ). Recent gravita- 
tional lensing studies ha v e discovered dar k substructures 
l|Dalal fc Kochanekl [200^ : iMao et al.|[2004D . so it appears 
ACDM is consistent with observations and we must find 
what physical mechanisms play the largest role in dark- 
ening small galactic halos. This leads to the other path, 
which is to consid er the effects of baryonic physics, such 
as st ellar feedback (iDekel fc Silklll986l:jMac Low fc Ferrara 
1999[) and UV ioiiisatio n (|Efstathioul Il992l : iQuinn et al. 
19961 : iBuUock et al.ll200oh . which might render many satel- 
lites dark. 

The four primary mechanisms that can remove mass 
from halos are: 

• UV ionisation: luminous objects emit UV radiation 
that ionises hydrogen and sets a background temperature 
above the virial temperature of the subhalo. 

• Ram pressure stripping: as a satellite passes 
through the hot halo gas, the incident gas pressure becomes 
stronger than the gravitational force of the satellite and gas 
is thus removed. 

• Stellar feedback: stellar winds and supernovae inject 
sufficient energy into the interstellar medium (ISM) of the 
small galaxy so that some or all of the ISM is ejected. 

• Tidal stripping: as a satellite orbits close to a larger 
host galaxy, the tidal forces become sufficient to remove ma- 
terial. Unlike the other three mechanisms mentioned above 
this is the only one that can remove coUisionless matter, 
namely dark matter and stars, as well as gas from a sub- 
halo. 

Previous efforts have been made at examining these 
mechanisms in detail. Early efforts were analytical due to 
the large dynamic range necessary to properly simulate sub- 
structures, but recent simulation s have allowed a closer look 
at satellites. iDekel fc Wool l|2003l ) compared careful observa- 



tions of many dwarfs with an analytical model based on the 
effect of supernova feedback and found that supernova feed- 
back defines the line between low and high luminosity dwarf 
galaxies. .Kravt sov et al.. (2004) used high resolution, cosmo- 
logical simulations to study the role of tidal stripping in the 
mass evolution of satellites. They concluded that the com- 
bined effect of tides and ionisation coul d produce a Milky 
Way-like satellite luminosity function. I Read et al.l (|2006t ) 
considered both supernovae-driven winds and ionisation in 
cosmological simulations and found that ionisation was criti- 
cal to ma ke their simulations agre e with observed luminosity 
functions. iGovernato et al.) (|2007) found in another series of 
cosmological simulations that UV background dramatically 
reduced the number of luminous subhalos, but that stel- 
lar feedback was required to make the simulated luminosity 
functions the same as those observed. More recent cosmo- 
logical simulations by Okamoto et al. (2010) have varied the 
strength of a kinetic supernova wind feedback to determine 
exactly how much energy is required to produce th e ob- 
served luminosity function. iKlimentowski et al.l (|2010l ) saw 
how tidal stripping determines a subhalo's baryori_ content 
and final morphology. Wadepuhl fc Springcl (20111) intro- 
duced black holes into their simulations and found that the 
black holes are not massive enough in subhalos to have an 
effect on their luminosities. They did, however, find that 
wind-driven galactic outfiows can reduce the number of high 
mass satellites, and cosmic rays can suppress the luminosity 
of low mass subhalos. Each of these models successfully fit 
the data by studying in detail one or two mechanisms, while 
we will consider all four within SPH simulations. 

Recent semi-analytic models also show some success at 
reproducing the observed satellite lumi nosity function. In 
these, only the more m a ssive subhalos (lOkamoto fc FrenkI 
(|2009l ), lGuo et all (|2010l ). lMacci6 et all (|2009 l)') retain stars. 
One semi-analyt ic model of an N-body simulation of a Milky 
Way-like halo (iLi et al.l bOOgl ) reveals luminous subhalos 
whose mass in dark matter spanned one order of magnitude, 
while the luminosity ranged over five orders of magnitude, 
matching observations. There were also many more dark 
matter-only subhalos present, whose mass spanned three or- 
ders ofnmgnitude^ 

iMaver et al.l (|2006l ) pointed out the importance of the 
combined effects of each mechanism. They simulated indi- 
vidual satellites falling into a static gravitational potential 
filled with hot, dense gas. In these simulations, tidal forces 
excite star formation and thus stellar feedback, as well as 
reshaping the gas distributio n so that it c a n be more easily 
stripped due to ram pressure. iMaver et all (120061 ) called this 
combination of processes "tidal stirring" and found that it 
can remove enough gas from dwarf irregulars to turn them 
into gas-free dwarf spheroidals. 

The purpose here is to discover which mechanisms tore 
the baryons off the subhalos, focusing on UV background, 
tidal stripping, ram pressure stripping, and stellar feedback. 
For the first time, we will explicitly track the causes behind 
the departure of individual gas particles from their subhalo 
in order to construct a comprehensive picture showing the 
relative strength of each gas loss mechanism. We only anal- 
yse the satellites inside the virial radius, so we expect the 
satellites we are analysing to be similar to dwarf spheroidal 
(dSph) galaxies, a population that dominates the satellite 
population of the Milky Way within Vvir- Dwarf spheroidals 
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are gas-poor (e.g. Fig 3 in lGrebel et"ai]|2003l ) containing as 
little as IO^Mq to undetectable amoun ts, but generally they 
continued forming stars until recently l|Skillmanll2007l ) . 

^establishes the background behind the tools we used 
in this work. iJ3] introduces the subhalo population of our 
host galaxy gl5784 at z = 0, while ijl] details the history of 
these subhalos and how we determine the causes of baryon 
loss. The concluding iJS] discusses the implications of our 
findings and areas for future work. 



2 METHOD 

We analyse the evolution of the subhalos of two galax- 
ies (gl5784 and g5664) from the McMaster Unbiased 
Galaxy Simulations (MUGS). The purpose of MUGS is 
to provide a sample of M* galaxies simulated using SPH 
at high resolution. A f ull de scription of MUGS can be 
found in IStinson et all (|2010l ). but we briefly summarize 
it here. The MUGS sample is chosen from a 50 Mpc 
volume of a WMAP3 ACDM universe {Ho=73 km s"^ 
Mpc~\ n^=0.24, nA=0.76, ntaryon^O.04, and (J8=0.79) 
ISpergel et al.ll2007l ). It consists of a random selection from 
the galaxies with halo masses between « 5 x 10^^ Mq and 
« 2 X lO^'^ Mq that did not evolve near to structures more 
massive than 5.0x IO^'^Mq within 2.7 Mpc. While this would 
have eliminated the Milky Way from our sample, there is no 
evidence for a past interaction between the Milky Way and 
M31 and so the Milky Way's satellite population should be 
unaffected by its near neighbor. The sample is unbiased with 
regards to angular momentum, merger history, and less mas- 
sive neighbors and it is hoped that the sample will reproduce 
the observed spread in galaxy properties. 

The only bias is random. The selected galaxies are sim- 
ulated with the commonly-used zoom technique that fo- 
cuses resolution on individual galaxies while maintaining the 
large scale torques necessary to give galaxies their angular 
momentum. The initial dark matter, gas and star particle 
masses are 1.1 x 10® Mq, 2.2 x lO'^ Mq and 6.3 x 10* Mq 
respectively. Each type of particle uses a constant gravita- 
tional softening length of 310 pc. Most of this paper will be 
dedicated to the subhalos of one of these galaxies, gl5784, 
which has a mass of 1.43 x 10^^ M!q. Its disk has a mass of 
3.27 x 10^° Mq and the bulge a n iass of 5.49 x 10^° Mr:^, 
based on kinematic decomposition (|Stinson et al.|[2010l l. We 
also utilize the 5.2 x lO^'^ Mq galaxy g5664 to see how the 
luminosity function changes depending on the presence of 
stellar feedback and UV background. 

Outputs were at most 214 Myr apart, especially at lower 
redshift, with irregular outputs at key times. Outputs were 
much closer together at high redshifts, typically 107 Myr 
apart. 

MUGS was run using the SPH code GASOLINE 
l|Wadslev et al.l 120041 ). gasoline includes low tem pera- 
ture metal cooling (described in IShen et al.l (|2010l ) and 
briefly here), UV background radiation, star formation, and 
physically-motivated stellar feedback. The metal cooling 
grid is con structed using CLOU DY (version 07.02, last de- 
scribed bv lFerland et al.l (|l998l ')'). assuming ionisation equi- 
librium. A uniform ultraviolet ionising ba ckground, adopted 
from Haardt & Madau (in preperation; see lHaardt &: Madaul 



(|l996l )). is used in order to calculate the metal cooling rates 
self-consistently. It starts to have an effect at z — 9.9. 



2.1 Star Formation and Feedback 

The star formation and feedba ck recipes are the "bl astwave 
model" described in detail in IStinson et all l|2006D . They 
are summarized as follows. Gas particles must be dense 
(nmin = 1.0 cm~^) and cool (Tmax = 15,000 K) to form 
stars. A subset of the particles that pass these criteria are 
randomly selected to form stars based on the commonly used 
star formation equation. 



dt 



^dyn 



(1) 



where M* is mass of stars created, c* is a constant star 
formation efficiency factor, Algas is the mass of the gas par- 
ticle spawning the star, dt is how often star formation is 
calculated (1 Myr in all of the simulations described in this 
paper) and tdyn is the gas dynamical time. The constant pa- 
rameter, c* , is tuned to .05 so that the simu lated isolated 
Model Milky Wa y used in lStinson et al.1 (|2006l ) matches the 
lKennicufr3 (|l998h Schmidt Law, and then c* is left fixed for 
all subsequent applications of the code. 

At the resolution of these simulations, each star parti- 
cle represents a large group of stars (6.32 xlO* Mq). Thus, 
each particle has its stars partitioned i nto mass bins based o n 
the initial mass function presented in iKroupa et all l|l993h . 
T hese masses ar e corre lated to stellar lifetimes as described 
in iRaiteri et al.l l|l996l ). Stars larger than 8 Mq explode 
as supernovae during the timestep that overlaps their stel- 
lar lifetime after their birth time. The explosion of these 
stars is t reated using the analytic m odel for blastwaves pre- 
sented in McKee fc Ostrikerl (Il977l ) as described in detail in 
IStinson et al. 1 20061) . While the blast radius is calculated 
using the full energy output of the supernova, less than 
half of that energy is transferred to the surrounding ISM, 
EsN = 4 X 10^° ergs. The rest of the supernova energy is 
assumed to be radiated away. 

To capture the behavior of clustered star formation, we 
stochastically determine when a star particle releases feed- 
back energy, 
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(2) 
(3) 



where Nsnii is the number of supernovae calculated to ex- 
plode during that star formation timestep, Nsnq is the 
"supernova quantum", i.e. the number of supernova re- 
quired per explosion (fixed at 30, the number of super- 
novae expected from the star particles in our simulation), 
and Nesn is the total number of supernova explosions that 
will have their energy distributed during a the star forma- 
tion timestep. If the probability, p, is greater than a random 
number, r, selected between and 1, Nsnq supernovae's 
worth of energy is released. This causes SN energy to be re- 
leased in quantized packets over the 35 Myr until the largest 
star remaining is < 8 Mq. 
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2.2 Group Finding: Amiga Halo Finder 



3.1 The Main Halo: gl5784 



In order to identify the host galaxy and its subhalos, we 
used the Amiga Halo Finder (AHF) (|Knollmann fc Knebd 
l2009h . AHF is based on the spherical overdensity method for 
finding halos. It is able to identify density peaks using an 
adaptive mesh algorithm. Once the density peaks are iden- 
tified, AHF cuts out halos (and subhalos) using isodensity 
contours. Particles belonging to subhalos are distinguished 
from those of the background halo using a simple unbinding 
procedure to determine whether the particles are gravita- 
tionally bound to the subhalo. We base our analysis on a 
minimum group size of 50 particles, which is 2.2 x 10^ Mq 
when the group only contains dark matter but could be less 
massive if it also contains gas and star particles. Our anal- 
ysis is restricted to those satellites identified by AHF as 
lying inside the virial radius (r^ir) of the halo. For gl5784, 
r^ir = 240 kpc, while for g5664, r^ir = 152 kpc. 



2.3 Merger Trees 

We traced the histories of each subhalo in the galaxy. First, 
we identified groups at every output 100 Myr apart with 
AHF and then traced the particles present in the subhalos 
at jz = back through the simulation including any gas 
out of which stars formed. For the sake of clarity, let us 
call the subhalo of interest "Alpha". At each output, we 
note every group that contains Alpha's particles. The group 
that had the largest number of Alpha's particles is set as 
Alpha's progenitor at that output. In this way, we trace the 
properties of each subhalo through time, including mass, 
distance from the host galaxy, and temperature. 

Because this method only depends on the number of 
particles, a subhalo can "jump" in position space between 
outputs, switching between subhalos with different groups of 
particles. Out of the higher mass subhalos, such behaviour 
was only observed in three of them. This does not affect our 
results because the jumps only occurred at high redshifts 
and between subhalos of comparable mass and position that 
were soon to merge. The pivotal point of analysis in each 
subhalo's history is its maximum mass, as will be show, in SJJI 
and it occurs well after any jumping behavior we see. Also, 
we will only look at the last time that a particle leaves its 
subhalo in order to prevent a single gas particle from being 
double-counted. Therefore, any event of gas being "lost" due 
to subhalo jumping will automatically be removed from the 
analysis. 



3 SIMULATED LUMINOSITY FUNCTIONS 

The first analysis we undertook was to compare the satel- 
lite luminosity function of gl5784 to observations. We found 
that the cumulative number is similar to the Milky Way's, 
though there was an excess of high luminosity satellites and 
so the shape of the cumulative functions did not match. We 
also re-simulated a smaller galaxy, g5664, with and without 
the UV background and stellar feedback to compare the ef- 
fects the presence these two mechanisms have on the subhalo 
population as a whole. 



AHF found 107 satellites inside r^ir of gl5784. In order to 
determine the luminosity of each subhalo, we treated each 
star particle as its own stellar population, where MUGS al- 
lowed us to model a population of stars with a distribution 
of ages. We based the brightness of the stars on the lu- 
minosity grid provided by CMD 2.1 (|Leitherer et al.lll999l : 
[Maxigo_et^al. 2008). Using the grid, we performed a bilinear 
interpolation over the stellar ages and metallicities of each 
star particle and then summed the luminosities of all the 
star particles in each satellite to derive a stellar magnitude 
for the satellites. We neglect the effects of dust extinction 
since dwarf ga laxies are low metallicity and rarely a ppear 
dust obscured (|Matedll998l : lLisenfeld fc Ferrarall 19981 ). This 
MUGS galaxy has an effective resolution of 2048^, and we 
will also compare gl5784 to a lower resolution run with an 
effective resolution of 1024^. 

Figure [1] shows the cumulative luminosity function of 
the host galaxy's subhalo population in the V-band at z = 0. 
The results show that down to My ~ —6, the dimmest 
subhalo in our simulation, the number of satellites is 23 
compared with 20 to 21 for the Milky Way at a simi- 
lar magnitude. Thus, when all the relevant baryonic pro- 
cesses are incl uded the order-of-magnitude mi s sing satel- 
lites problem (|Moore et al.l 1 19991 : iKlvpin et al] Il999l ') dis- 
appears as has b een seen in other simulations of compa- 
rable resolution ([Okamoto et all l2010l : iKnebe et all I2OIOI : 
IWadepuhl fc SpringellbonF T 

For comparison to observations we include two lines for 
reference. This first is a recent compilation of the classi- 
cal satellites and the new ultra-faint dwarf galaxies from 
ToU erud et al] (|2008l l. The census of ultra-faint dwarfs is 
certainly incomplete, both due to the limited sky coverage 
of the Sloan Digital Sky Survey and the difficulty in detect- 
ing very faint galaxies. These galaxies extend to lower lumi- 
nosities than our simulations can resolve, where our faintest 
satellite has one star particle. The resolution of our simula- 
tions is not sufficient to make robust predictions regarding 
the new classes of ultra-faint dwarfs. 

The second is iKoposov et"al] (|2008l ) who modelled these 
effects and derived a corrected luminosity function, — 
10 X io° i(*^v+5) fgj. _;^g ^ j^j^ ^ „2, that would rep- 
resent a theoretically complete set of satellite galaxies; we 
have also plotted this function in Figure [T] Our host galaxy 
gl5784 has a mass close to the Milky Way's, so it is rea- 
sonable that its cumulative number of luminous subhalos 
is comparable to the Milky Way's, as indeed it is. In other 
words, our simulated galaxy does not suffer from the missing 
satellite problem. The major difference between the simula- 
tions and observations is an excess of brighter satellites in 
the simulations, and that our star formation recipe may form 
too many stars. This causes an extra "knee" in the shape of 
our luminosity function. 

The dash-dotted line in Figure [1] shows the luminosity 
function of a lower resolution simulation exactly the same 
as gl5784, but with half the spatial resolution and an initial 
gas particle mass of « 10® -'^0- From this it is clear that 
decreasing resolution decreases the number and luminosity 
of subhalos. The lowest luminosity subhalo, at Mv ~ —8.2, 
contains a single star; in our higher resolution run this mag- 
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Figure 1. The cumulative V-band luminosity function of the host 
galaxy's subhalos at z = 0, in s olid black. Ob servational data of 
the Milky Way as gathered by iTollerud et al . ( 2008 ) is shown 
in da shed red, while a theo retical completion-corrected function 
from IKoposov et alj ||2009| ') is shown in dashed blue. The black 
dash-dotted line is the luminosity function from a lower resolution 
run of the same galaxy where the initial gas particle mass is 10® 

Mq. 

(Note that the observed ultra-faint dwarfs extend to a much 
lower luminosity). Within the resolution limit, all three 
functions lie within an order magnitude of each other. 



nitude corresponds to 10 star particles. We discuss resolution 
effects in ^3.31 

Figure O shows the baryonic mass of each subhalo as 
a function of total mass at 2: = 0. The dashed line shows 
where subhalos that obey the cosmic baryon fraction would 
lie. The subhalos contain systematically fewer baryons below 
5 X 10* M(7^. This falloff is simila r to the low mass dropoff 
found bv lMcGaugh fc Wolj (|2010l ) in the observed baryonic 
TuUy-Fisher relationship. In these lower mass halos, baryons 
are preferentially stripped. As we shall see, these halos have 
also lost a significant amount of dark mass. However, almost 
all the lower mass halos also have fewer than 10 baryons. 

Additionally, below 2 x 10® Mq there are many halos 
that contain no baryon particles according to our simula- 
tion's resolution at all and are thus dark. Of the 23 satellites 
that contain baryons, only 10 contain gas with the maxi- 
mum gas fraction being 4 % of the total mass. This fraction 
might seem low if compared to high gas fractions found in 
isolated dwarf irregular galaxies (|Geha et al.ll200^ ), but the 
subhalos considered here are all within r„ir and are more ap- 
pr opriately compared with the dwarf spheroidals presented 
m iGrebel et al] (|2003D . which contain little gas. We must 
stress that when we say that a subhalo "has no baryons" we 
mean in the sense of our simulation's resolution there are no 
baryon particles. Actual dark subhalos will always have at 
least some trace quantities of gas. 

Figure [3] shows the cumulative mass function at 2: = 
for all the satellites as well as the subset that formed stars. 
The total satellite mass function is simil ar to th e coUision- 
less, dark ma tt er-on ly simulations of .Moore et al.l |l999) and 
iKlvpin et all (|l999l ) while the mass function of luminous 
satellites is closer to what is observed. There is about an 
order of magnitude more subhalos in total than those that 
contain baryons. How these dark subhalos lost their gas and 
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Figure 2. Baryon mass versus total mass for our subhalo popu- 
lation. The symbols correspond to mass to light ratio (total mass 
divided by baryonic mass), and the dashed black line is where 
subhalos with the cosmic baryon fraction would lie. The lower 
bound of the mass-to-light ratio is inclusive. The horizontal dot- 
ted line is the most luminous of the satellites with fewer than 
10 baryon particles, corresponding to our resolution limit. Most 
striking here is that the luminous subhalos, the brightest of which 
follow the power law, and dark subhalos overlap in total mass. 
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Figure 3. The cumulative mass function of our host galaxy's 
subhalos at z = 0, divided into all subhalos in blue, and those with 
baryons in red. The total number of subhalos is over a hundred, 
but the total number containing baryons is almost an order of 
magnitude fewer. 



their stars, if they ever had any, is the key to understanding 
the missing satellites problem. We will investigate gas re- 
moval mechanisms, the reason some subhalos are light and 
some are dark, and why the only 23 subhalos that are lumi- 
nous are not also the 23 most massive in ^ 



3.2 Effects of Stellar Feedback and Ultraviolet 
Background: g5664 

One test to delve deeper into how the cumulative mass func- 
tion compares with dark matter-only simulations while the 
cumulative number of luminous satellites is consistent with 
observations is to re-simulate a galaxy with and without the 
UV and stellar feedback. We re-simulated a second, smaller 
host galaxy, g5664, with a mass of 5.2 X 10" M0, three 
times with different baryonic physics: 
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(a) with the standard MUGS simulation including UV 
and stellar feedback, 

(b) with UV but no stellar feedback, and 

(c) with neither UV nor stellar feedback. 

This galaxy contained half the number of gas and star par- 
ticles as gl5784 and thus was faster to run, which is why 
it was chosen for the parameter comparison. Since gl5784 
is more massive and has a larger r^ir, it will contain more 
and likely more massive satellites. Since the UV background 
and stellar feedback more strongly effect low mass galax- 
ies, gl5784 will contain more luminous satellites. However, 
this should not affect the relative comparison of the differ- 
ent resimulations of g5664. The relationship between these 
mechanisms and the mass of the satellites will be explored in 
^ in addition to how these mechanisms affect the satellites 
throughout their history. 

Figure U shows the cumulative luminosity function for 
the three simulations at z = 0. ft confirms the reduction 
in number of satellites. There about one-fourth as many in 
g5664 as in gf 5784. Nearly every satellite in the simulation 
run without feedback and UV (c) contains stars. Conversely, 
many fewer satellites contain stars in the simulation that in- 
cludes feedback (a). It is apparent that the UV ionisation 
plays a large role in stopping stars from forming in many 
subhalos. When stellar feedback is added, there is little ef- 
fect on subhalos brighter than My ~ —15, but fainter sub- 
halos are only populated with a few stars whose feedback 
was effective at eliminating star formation for the rest of 
the simulation. The My « — 15 threshold is similar to the 
flattening seen in Figure [T] for both the simulated gl5784 
and the observed Local Group mass function. 

The UV-only simulation (b) contains several medium 
luminosity subhalos but none that are extremely faint, stop- 
ping at Mv ~ —12. It is curious then that simulation 
(a), with both feedback and UV, lacks satellites between 
—8 < Mv < —15, but contains two very faint satellites 
at Mv ~ —6. The two luminosity functions diverge at 
Mv ~ —15, which corresponds to about 10^ star particles 
in all three runs of g5664. The origin of this situation is un- 
clear from inspection of the luminosity functions alone. In 
SjH however, when we track the mechanisms of baryon loss 
through time an explanation for this phenomenon will arise 
showing that stellar feedback preferentially strips medium 
mass satellites. 



3.3 Resolution 

Satellite galaxies are the most difficult objects to resolve 
in simulations, and so they show the strongest resolution 
effects. Figure [T] begins to show how resolution can affect 
satellites. It is harder to detect low mass satellites in lower 
resolution simulations, and effects such as tides and ram 
pressure stripping are more accurately captured as resolu- 
tion increases. 

Other authors have discussed the resolution effects on 
their satellite luminos ity functions in simulations similar to 
ours. iLibeskind et al ] ((2007) ran simulations at a lower res- 
olution than ours, studying the satellite systems of several 
ACDM galaxies with their gas resolved to « lO'' Mq . They 
compared their luminosity functions to semi-analytic mod- 
els of resolution over four different orders of magnitude and 
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Figure 4. Cumulative luminosity function for g5664, comparing 
runs with three different conditions on the baryonic physics. The 
run with UV background and stellar feedback is in red (a). The 
run with U V background, but no feedback is in green (b) . The run 
with neither UV background nor stellar feedback is in blue (c). 
Both feedback and the UV background reduce the total number 
of luminous subhalos. 

found a convergence at Mv ~ —12. lOkamoto et al.l (|2010l ) 
studied the effects of several feedback models on the satellite 
populations on host galaxies around the same mass as our 
gl5784 and with similar gas particle masses. Matching the 
circular velocity of their satellites to a power law, they con- 
cluded that their satellites were well res olved down to satel- 
l ites w ith at least 10 star particles. Wa depuhl fc Springell 
(j201ll i compared a low resolution simulation (whose gas par- 
ticle mass is « 10^ Mq corresponding to our high resolution 
simulation) to their high resolution simulation with gas mass 
of ~ 10 Mq and found them convergent up to Mv ~ —8. 

We compared our high resolution simulation of gl5784, 
with initial gas particle mass « 10^ Mq, to a low resolution 
simulation, with initial gas particle mass « 10*^ Mq, and 
found that their luminosity functions were similar in their 
region of overlap down to Mv ~ —8.2. The major difference 
is that at higher resolution there appears to be a slight over- 
abund ance of highly luminous subhalos. IChristensen et al] 
(|2010l ) examined the effects of resolution on galaxies rang- 
ing from 10^^ to 10^ Mq and found that 10* gas particles 
at each galaxy's maximum mass are required before the star 
formation recipe used here converges. 

Subhalos that formed stars in our high resolution 
gl5784 simulation had their star formation drop off at about 
10* star particles at z = 0, corresponding to the number of 
gas particles being at least ~ 5 x 10'^ at maximum mass. 
The star count corresponds t o Mv ~ —15 in our high res- 
olution simulation. Based on IChristensen et al.l l|2010l ). we 
expect to see a decrease in star formation across this range. 
However, it is difficult to conclude that the reduction in star 
formation is solely due to resolution since it corresponds to 
the mass range between 10^ and 5 x 10^" Mq where feed- 
back has its strongest effects. Our low resolution simulation 
of gl5784 only has one subhalo at Mv ~ -17 with ^7x10^ 
star particles, and only two more around Mv ~ — 16 with 
« 10^ star particles. The discrepancy between high and low 
resolutions is explained by the lack of particles in the low 
resolution run. 

Some of our subhalos contain many fewer particles 
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than IChristensen et al] l|2010l ) suggest is necessary to re- 
solve star formation, but we include them for completeness. 
We will still show and analyze all luminous and dark sub- 
halos with comprehensive h istories (as explained ea rly in 
311 . With higher resolution, IChristensen et al.l (|2010l ') sug- 
gest that these satellites would contain slightly more stars 
since gas could form a disk and become denser than the 
star formation threshold. However, we note that it is hard 
to draw conclusions without the higher resolution simula- 
tions because of the non-linear interaction of star forma- 
tion and stellar feedback, particularly considering that the 
present simulation over-predicts the number of highly lumi- 
nous satellites (Stinson et al. 2010). 



4 GAS LOSS MECHANISMS 

The previous section showed that both stellar and UV feed- 
back play a significant role in the evolution of satellite gala^c- 
ies. In order to gain an overall picture of how stellar and UV 
feedback work with ram pressure and tidal stripping, we will 
focus on gl5784 and study the detailed, particle-by-particle 
mass loss history of its subhalos. 

We traced the evolution of 85 of the 107 subhalos identi- 
fied by AHF in gl5784. We were unable to perform a detailed 
trace of every subhalo for two reasons: 

(1) 17 low mass subhalos did not maintain 50 member 
particles throughout the simulation, which is 2.2 X IO'^Mq 
when the group only contains dark matter but could be less 
massive if it also contains gas and star particles. 

(2) 5 halos were spatially coincident with another subhalo 
during one output, and thus appeared to gain a large amount 
of mass. Since our analysis focused on cumulative baryon 
loss and the subhalos' maximum mass these sudden spikes 
in mass would invalidate any results including those halos. 

10 of the 85 subhalos we analysed contained both gas 
and stars at z = 0, with total masses ranging from 5.6 x 
10** Mq to 3.3 X 10^° Mq. Of the subhalos that did not 
have gas at 2 = 0, 17 formed stars at some point but only 
13 of these retained them until z = 0. This leaves a total of 
23 luminous subhalos at z = 0, of which 13 have more than 
50 baryon particles. 

Figure [5] shows examples of the time evolution of these 
subhalos' mass and distance to the host at each output, 
representing an upper limit on the subhalos' closest distance 
to the host. The top panel shows subhalo (a) that retains 
gas and stars at 2; = 0, the middle panel shows subhalo (b) 
that retains stars but not gas at 2 = 0, and the bottom 
panel (c) shows a dark satellite that has no baryon particles 
at 2 = 0. Close passages to the center of the main halo most 
strongly removes gas and to a lesser extent, dark matter 
and stars through tidal stripping. Stars tend to sit at the 
centre of the subhalo and are less vulnerable to stripping 
than the dark matter around the subhalo's exterior. That 
is, subhalos that form stars before they lose their gas retain 
those stars until 2 = 0. If stars are not formed before their 
first close passage then the subhalo will never form stars 
and becomes a dark satellite. The quantity of gas lost is 
also determined by the subhalo's proximity to the host. For 
example the luminous first and second subhalos (a) and (b) 
start out with similar mass, but the latter has a pericentre 
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Figure 5. The time evolution of the distance to the host and 
the mass of different components of three subhalos. The black 
line corresponds to the distance to the host at each output on the 
right-hand axis on each plot. The coloured lines correspond to the 
masses on the left-hand axis: dark matter (blue), gas (red) and 
stars (green). The top subhalo (a) is massive enough to retain 
both gas and stars to ^ = (final mass: 3.9 X 10^ Mq). The 
centre subhalo (b) retains stars but not gas at 2 = (final mass: 
2.0 X 10* Mq). The bottom subhalo (c) loses all its baryons before 
2 = 0, making it a dark satellite (final mass: 7.2 X 10® Mq) despite 
the fact that it ends up heavier than the luminous subhalo in the 
centre panel. 



that is more than twice as close to the host and as a result, 
it loses all its gas by 2 = 0. Central location is not the full 
story. The dark satellite (c) is farther from the host galaxy 
than either of the luminous ones and has a higher mass at 
2 = than the subhalo (b), but it never forms stars and 
its mass in gas is much lower before being lost early in the 
simulation. There are more factors than a subhalo's mass at 
2 = and its proximity to the host galaxy that decide if it 
is luminous at 2 = 0, which we will explore in this section. 
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4.1 Ultraviolet Background 

Quantifying iiow mucli gas a sublialo loses due to the UV 
background involves some degree of subjectivity. Instead of 
merely counting gas that was in the subhalo, it will involve 
determining how much gas was never in the subhalo but 
ought to have been. 

We confine our analysis of the evolution of subhalos to 
the time after z = 10.8, when the first halos get larger than 
the minimum group finder particle limit and are virialized. 
One of the earliest mechanisms that removes gas from sub- 
halos is the reionising UV background radiation that is emit- 
ted from the first luminous objects. This radiation makes its 
first impact at z = 9.9 in these simulations. 

Most of the subhalos cannot be identified until a few 
outputs after z = 9.9. To enable the analysis of the effects 
of UV radiation, every dark matter particle within r^ir is 
matched with a "twin" gas particle at the initial condi- 
tions. In our halo-by-halo analysis, dark matter twin par- 
ticles are defined as subhalo members at the time their sub- 
halo reaches its maximum mass. The evolution of the gas 
twins of these member particles is then traced from the ear- 
liest output onwards. The ensemble of the twin gas particles 
are referred to as the "background gas" later in this section. 

The mass lost due to reionisation is defined as the gas 
that had a dark matter twin in a given subhalo, but it- 
self was never contained in that subhalo. The reason that 
this gas was not in the subhalo is that subhalos are unable 
to contain gas with a temperature higher than the subhalo 
virial temperature (T^ir)- Low mass subhalos will thus con- 
tain dark matter without its gas twin. Note that the twins 
are calculated from dark matter in the subhalo at the time 
of the subhalo's maximum mass and not reionisation. This is 
important because the way subhalo tracing works, only one 
subhalo is labelled as "the" subhalo at any given timestep. 
However, the larger subhalos reach their maximum mass as 
a result of the merger of several smaller subhalos. Determin- 
ing the UV loss at a subhalo's maximum mass accounts for 
loss in each of the individual subhalos. Figure |S] shows that 
satellites reach their maximum mass prior to being captured 
by the main halo, typically immediately before. 

Subhalos that did not virialize before reionisation gen- 
erally do not contain gas. The subhalos that virialize after 
reionisation, but contain gas that goes on to form stars, do 
typically have gas without a dark matter twin in the sub- 
halo. In the face of such a mismatch between gas and dark 
matter particle members we developed the following analy- 
sis to justify why the amount of background gas that ends 
up in the subhalo correlates to UV loss. 

In order for a gas particle to have enough energy to 
overcome the potential of the subhalo and escape, its tem- 
perature must be greater than Tvir, where; 

2GjJ,mpM subhalo 



3kRs 



(4) 



where G is the gravitational constant, /x is the mean molecu- 
lar weight (where the typical value is 0.6 for ionised gas), mp 
is the proton mass, k is Boltzmann's constant, and M subhalo 
is the subhalo's mass. Rsubhaio is the distance between the 
halo's centre of mass and the most distant member particle. 

Figure [7] shows the evolution of the background mean 
temperature (twin particles) and Tvir in halos (b) and (c) 
from Figure [5l The subhalo in the top panel (b) retains its 
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Figure 6. For all our subhalos, the time of maximum mass ver- 
sus the time they join the main halo. Most of them reach their 
maximum mass shortly before becoming substructure. The dotted 
line marks where the time of maximum mass is equal to the sub- 
structure time, showing that none of the subhalos achieve their 
maximum mass after becoming substructure. 



baryons through reionisation because of its higher mass and 
consequently higher Tvir- While halo (c) is more massive at 
z = than (b), it is dark because the twin gas of its dark 
matter was hotter than Tvir - The subhalos that form earliest 
capture the most gas and form the most stars. 

Figure |8] shows how the ratio of T^ir with the mean gas 
twin temperature effects the final baryon fraction. Specifi- 
cally, we compare the maximum ratio of Tvir to background 
gas mean temperature versus the baryon fraction at two 
times: z = and the time at which this temperature ratio 
is a maximum. The subhalos with T^ir /Tbackground > 1 at 
some point are more likely to be near the cosmic baryon frac- 
tion. With one exception, the subhalos that contain baryons 
at 2 = are a subset of these. The one exceptional satellite 
started out near to the host and gathered enough gas early 
in its formation to form enough stars to remain luminous 
over its subsequent several close passages to the host, even 
though it lost all its gas on its second passage. Since none of 
the subhalos that never contained their twins held onto any 
of their baryons, it appears that the "twin" particle analy- 
sis is robust. Therefore, we can safely define gas lost due to 
UV background as the mass of twin gas particles that never 
entered subhalo. 



4.2 Ram Pressure Stripping 

One of the difficult aspects of this study is that the mass 
loss mechanisms involve the interaction between two gas 
phases with significantly different propert ies. This ci rcum - 
stance is one where SPH struggles. Ag ertz et al] (|2007l ) 
showed that SPH has trouble modelling ram pressure strip- 
ping, particularly when the Kelvin-H elmholtz time is im- 
portant. However, iMaver et all (|2006l ) showed that SPH can 
model stripping when Tdyn < t kh, and our g l5784 satellites 
typically fall into this regime l|Maver et al.l i2006 show that 
T"KH ^ 4 Gyr when they reach their minimum at pericentre, 
compared to satellite dynamical times of Tdyn « 0.2-6 Gyr). 

With thes e caveats in mind, w e do a classical ram pres- 
sure analysis (|Gunn fc Gotdll972l ) to see how close these 
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Figure 8. The baryon fraction versus the maximum ratio be- 
tween the virial and background gas temperatures that a subhalo 
ever obtains over its lifetime. The horizontal dashed black line 
represents the cosmic mean, 0.17, while the vertical dashed black 
line separates the subhalos that achieved a T„i,. higher than the 
background gas temperature on the right from those that did not 
on the left. The asterisks represent subhalos that have ten or more 
baryon particles at 2 = 0, while Xs represent subhalos that have 
fewer than ten baryon particles at z = 0, corresponding to the res- 
olution limit found in Figure [T] Red shows the subhalos' baryon 
fraction at the time of their maximum ratio of T^ir /Tbackground^ 
while blue shows the subhalos' baryon fraction at z = 0. With 
one exception, every subhalo that is luminous at 2 = is to the 
right of the vertical line. 



Figure 7. Evolution of background gas mean temperature of the 
background twin gas (blue) and T^ir as given by Equation|4](red) 
for subhalos (b) and (c) from Figure [5] Temperature corresponds 
to the left axis, and the subhalos' orbits (black, solid) to the 
right. The vertical dashed line corresponds to z = 9.9 when the 
UV background is turned on. Even though the subhalo on the 
bottom is more massive at 2 = 0, it is unable to retain as much 
gas as the lower mass subhalo on the top due to the fact that its 
background gas is consistently hotter than T^ir- 



simulations come to reality. In order to quantitatively mea- 
sure the efTect ram press ure stripping has on our subhalos, 
we use the criterion from lGrebel et al] (|2003h : 



^ram ~ Phhg^subhalo ^ q 



(5) 



where Pram is the ram pressure, phhg is the gas density in 
the hot halo gas around the subhalo, Vsuthaio is the subhalo's 
velocity relative to the host galaxy, Usubhaio is the velocity 
dispersion of the gas in the subhalo, pgaa is the average den- 
sity of the gas in the subhalo. Here the subhalo's velocity 
dispersion is defined as: 

2 3 GMs^hhalo /„\ 

'^subhalo — v"-' 

^ ^subhalo 

where Mgubhaio is the subhalo's mass, and Rsubhaio is the 
subhalo's radius. 

The gas density of the hot halo gas around each halo is 
defined as the average density of the n nearest gas particles, 
where n is twice the number of particles in the subhalo to a 
maximum of 4000. The density and temperature structure 
of the gaseous halo of gl5784, through which the subhalos 
pass, is shown in Figure [9] with substructure removed. 

Every gas particle that leaves during the timesteps when 



0) 

D 

D 

QJ 
CL 

E 



10' 




10000. 



1000 



K) 

o 

Q. 



c 



Figure 9. The gas temperature (red) and density (blue) profiles 
of gl5784 at 2 = for the outer disk at r ^ 30kpc. Power laws 
fit to p oc r~^''^ (dashed purple) and T oc r-— 0-73 (dashed orange) 
are shown for reference. 



the ram pressure exceeds the internal pressure of the halo 
is classified as leaving due to ram pressure stripping except 
for those particles that quahfied for stellar feedback ( ^4.4|l . 

However, we also generated movies of twelve low-mass 
subhalos and visually inspected them to see if ram pressure 
or tides removed their gas. Figure [10] gives an example of a 
subhalo that loses its gas due to ram pressure. As a subhalo 
approaches the host galaxy, it enters the hot halo gas. In 
low enough mass subhalos, its gas gets left behind, cleanly 
separating from its dark matter. Figure [10] shows that the 
gas maintains the shape of the subhalo for a several tens of 
Myrs. This contrasts with the signature of tidal stripping 
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Figure 10. How ram pressure stripping removes gas from a 
1.7 X lO'^ Mq subhalo between z = 0.26 (left) and z = 0.16 
(right). The dark matter (dark green), gas (hght green) and stars 
(turquoise, none present here) are marked at the time of max- 
imum mass, while the bottom layer (brown, covered the other 
layers) is the subhalo at the present output. Ram pressure strip- 
ping has affected this subhalo, removing all of the gas it pocessed 
at the time of maximum mass, but without strengthening the 
ram pressure stripping diagnosis by a factor of 10 the individual 
gas particles would not be marked as ram pressure stripped. The 
background colours denote the temperature of the surrounding 
gas, ranging from dark blue (colder, 10^'^ K) through to white 
(hotter, 10^ K). 





Figure 11. Pressure versus time for subhalos (b) and (c) from 
Figure [5] Internal pressure is given by the solid red line, ram 
pressure from the hot halo gas is the blue dashed line, and ram 
pressure strengthened by a factor of 10 is the green dashed line. 
Visual inspection of the subhalos revealed that the one in the left 
panel was affected both by ram pressure stripping and later tidal 
stripping, while the subhalo on the right was affected predomi- 
nantly by ram pressure stripping. 



where tidal forces elongate the matter ahead and behind 
the satellite's orbit. Typically, ram pressure stripping occurs 
farther out from the host than tidal stripping. 

Our definition for ram pressure underestimated the ef- 
fect that is visible in the simulations. This seems to be a 
numerical effect and requires a multiplication of the calcu- 
lated ram pressure by a factor of 10. The factor of 10 may 
compensate for the use of the mean satellite gas pressure 
instead of the pressure in the outer regions where particles 
are getting removed. Figure [TT1 shows how the subhalo's gas 
pressure and ram pressure evolves in two subhalos. The fac- 
tor of 10 increases the ram pressure so that it is comparable 
with the subhalo's gas pressure. In the future, a compari- 
son of higher-resolution simulations, as well as grid codes, 
to ours would be useful to see whether they exhibit results 
closer to the analytic determination. 



4.3 Tidal Stripping 

Next, we examine mass loss due to tidal forces using a simple 
spherical approximation for the host halo and the satellite. 
For each subhalo, the tidal force particles feel from the host 
is compared with the force they feel from their internal grav- 
ity. The subhalo's tidal radius is the place at which its self 
gravity is less than the tidal force of the host galaxy (e.g 
iHavashi et ahl (|2003l )'). Due to Newton's Theorem, assum- 
ing spherical symmetry, one need only consider the mass 
interior to a given particle. The gravitational force on the 
particle from inside the particle's orbit is 



F. 



subhalo 



ubhalo 



(7) 



where r is the distance from the particle to the centre 
of its subhalo, mparucie is the mass of the particle, and 
Msubhatoir) is the subhalo's mass interior to r. The tidal 
force that the particle feels from the host, if the host galaxy 
is approximated by a point mass and assuming all its mass is 
contained inside the satellite, is the differential pull between 
the particle's position in the satellite and the satellite's cen- 
tre: 



SFtidal 



-2GMhost'mparticler 
-"•host 



(8) 



where Rhost is the distance between the subhalo and the 
host, and Mhost is the mass of the host. Therefore, the con- 
dition for when the particle feels a greater tidal force than 
gravitational from its own subhalo is given by: 



Msubhalo{r) 2Mhoa 



^hos 



(9) 



A particle that passes this test qualifies for tidal stripping. 
We emphasize that this a spherical approximation, given 
that the subhalos occasionally have their shape distorted, 
though the distortion is usually symetrical. We did try vary- 
ing the strengh of the tidal force, as we had done for ram 
pressure stripping, but found that it did not change results 
much and so this method has some degree of robustness. 

Figure [121 shows an example of a subhalo being tidally 
stripped of its gas. The tidal force pulls material out in lead- 
ing and trailing arms. Additionally, material is stripped off 
the outside first before the material on the inside. 

Regarding the possibility of tidal stirring. Figure [5] 
shows that there was not significant star formation follow- 
ing close passages of satellites. During these close passages, 
dark matter is often stripped and a large fraction of gas 
is stripped. It is possible that these simulations are too 
low of a resolution to model tidal stirring like was seen in 
iMaver et al.. (,2006, ). 

4.4 Stellar Feedback 

Often, much credit for the removal of baryons is given to stel- 
lar feedback (iDekel &: "Siii3ll986l : [Mac Low fc Ferrara|[T999l : 
iDekel fc Wodl2003l ). Supernovae release large amounts of 
energy into the ISM, which can be sufficient to liberate gas 
from halo potential wells. Others have argued that the cou- 
pling between the stellar feedback and the ISM is insuffi- 
cient to remove a significant amount of gas from subhalos. 
Determining how much gas stellar feedback removed in the 
simulations proved to be a challenging task. 
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Figure 12. How tidal stripping removes gas from a 2.0 X 10 Mq 
subhalos between z = 1.0 (left) and z = 0.8 (right). The dark 
matter (dark green), gas (light green) and stars (turquoise) are 
marked at the time of maximum mass, while the bottom layer 
(brown, covered the other layers) is the subhalo at the present 
output. Note how the tidal stripping has already begun at z = 1.0, 
and by z = 0.8 the tidal tail is prominent. The dark matter has 
been compressed, while the stars sit more safely in the subhalo's 
centre. The background colours are gas temperature as in Figure 

One signature of stellar feedback in these simulations is 
that a gas particle has its cooling turned off, which allows 
it to maintain its high temperature due to the stellar en- 
ergy release. When it gains sufficient kinetic energy, it can 
escape the gravitational potential of the subhalo. We found 
this method of tracking stellar feedback to be very limiting, 
however. Outputs were limited to approximately one every 
200 to 100 Myr. Cooling is typically shut off for 50 Myrs, so 
between 1/2 and 3/4 of particles whose cooling was turned 
off would be missed by simply counting particles whose cool- 
ing was shut off during an output. 

Another signature of stellar feedback is the release of 
metals into the surrounding interstellar medium. In our sim- 
ulations, a star particle releases about 300 Mq of metals 
from type II supernovae to the nearest 32 gas particles. The 
ejection is smoothed so that gas closer to the star receives 
more metals than particles further away. A typical gas par- 
ticle will receive a few Mq in metals from a star particle. 
Gas can also receive metals from other gas particles through 
diffusion. These metal transfers are typically < 1 Mq. 

So, when gas had an increase in metals of ^ 5 Mq, it 
is likely that it was in the neighborhood of stellar feedback 
and we classify it as having been lost from the halo due to 
stellar feedback. 

This is a conservative estimate because there may also 
have been cases where gas directly heated by stellar feedback 
acquired sufficient pressure to push out different gas, a pro- 
cess called "mass-loading" . Thi s is common in dwarf galaxies 
(IDalla Vecchia fc Schavdl2008l ). Stellar feedback could have 
augmented another mechanism like ram pressure stripping 
and gotten unlabelled as such. 

4.5 A Combination of Mechanisms 

Using all the techniques described above, we now present a 
summary of which processes dominated accretion and loss 
of gas. 

One confounding effect happened when massive subha- 
los pass through the pericentre of their orbit. Subhalos tem- 
porarily accreted a small quantity of gas and quickly lose 
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Figure 13. Mechanisms for baryon loss over time for the massive 
subhalo (a) (maximum mass: 7.4 X 10^ Mq, final mass; 3.9 x 
10^ Mq) from Fi gure [5] Turquoise is the cumulative gas lost 
due to the UV background, red is the cumulative gas lost due 
to ram pressure stripping, dark green is the cumulative gas lost 
due to tidal stripping, light green the is the cumulative gas lost 
due to stellar feedback, purple is the cumulative gas lost due to 
undetermined causes, orange is the cumulative lost stars, blue is 
the gas at the current timestep, and yellow is the stars at the 
current timestep. 

it, possibly a numerical effect. Because of this, we did not 
count mass loss of particles that entered and left subhalos 
after they reached their maximum mass. 

Some particles left the subhalo more than once. In order 
to avoid double counting, only the last method by which a 
particle entered or left its subhalo was counted. We noted 
which gas particles were converted into stars, so that we 
could identify what portion of the gas mass decrement was 
due to star formation and what portion was due to gas leav- 
ing the subhalo. The loss of many particles could not be 
classified, so mass loss is often classified as "other", empha- 
sizing the schematic nature of this method. 

Figures [TH] to [TS] show the evolutionary history of the 
same subhalos as shown in Figure (5] These histories are 
coloured to indicate the amount of gas lost due to each of the 
mechanisms described above, as a fraction of each subhalo's 
maximum mass. Most of the baryons for the most massive 
subhalo (a) turn into stars quickly, and, other than the initial 
UV ionisation which prevents a significant amount of gas 
from being captured, any gas that is lost is usually due to 
tidal stripping or stellar feedback. After UV ionisation the 
medium mass subhalo (b) that retains its stars loses its gas 
largely due to stellar feedback with smaller contributions 
by tidal and ram pressure stripping, though there is a large 
number of unclassified particles as well. The medium mass 
subhalo (c) that ends up as a dark satellite lost almost all of 
its gas due to the UV ionisation, with the small remainder 
being stripped by ram pressure. 

Figure [16] combines all plots of type Figures [13] to [15] 
for all subhalos and shows as a fraction of each subhalo's 
maximum mass the gas lost due to each mechanism a,t z = 
0. These mass loss fractions are plotted as a function of 
maximum mass rather than final mass because the sequence 
of mechanisms appears more clearly (as shown in Figure [21). 
Figure [16] shows that the massive subhalo (a) in Figure [13] is 
no aberration. It is common for the most massive subhalos to 
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Figure 14. Mechanisms for baryon loss over time for the medium 
mass subhalo (b) (maximum mass: 4.2 X 10^ M0i final mass; 
2.0 X 10** Mq) from Figure [5] The colour scheme follows Figure 
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Figure 15. Mechanisms for baryon loss over time for the medium 
mass subhalo (c) (maximum mass: 1.4 X 10^ Mq, final mass: 
7.2 X 10* Mq) froni Figure [5] The colour scheme follows Figure 



efficiently form stars and for tidal stripping to play the most 
important role in their mass loss. Medium mass subhalos 
tend to be more dominated by stellar feedback, since they 
are massive enough to form stars but light enough that they 
are more susceptible to losing their gas. Following Figure 
1151 lower mass subhalos lose significant mass due to UV 
reionisation and then much of the remaining mass is stripped 
by ram pressure. 

A dichotomy of evolutionary scenarios appears in Figure 
1161 Halos less massive than ~ 2.0 x 10^ Mq lose their gas due 
mostly to UV reionisation. Higher mass halos formed stars, 
lost less mass to reionization, and lost gas due to a variety 
of other mechanisms. This distinction disappears when the 
satellites are classified by their final mass as in Figure [2] 
where the populations of baryonless and luminous subhalos 
overlap in terms of total mass at a = 0. 

Table [1] summarizes the results of Figure [16] by dividing 
the subhalos into three categories (ones with gas and stars 
at z = 0, ones without gas at z = but had stars, and 
ones that never formed stars and have no gas at z = 0) . UV 
ionisation is the most prominent for all subhalos. For the 




Max MasSp,<,|„(MQ) 

Figure 16. Mechanisms for baryon loss at z = Q for all subhalos, 
as a function of the maximum mass that a subhalo was able to 
achieve. The colour scheme follows Figure 1131 The various loss 
mechanisms are cumulative over time, while the total gas and 
stars are for the current time. All values are a fraction of the 
subhalo's maximum mass. The luminous subhalos show an excess 
of baryons over their lifetime, above the cosmic mean of 0.17. 
The cutoff between luminous and dark satellites happens at about 
2.0 X 10^ Mq. 



most massive subhalos tidal stripping followed, while stellar 
feedback and ram pressure slipping had little impact. For 
the subhalos massive enough to form stars at some point 
but which did not retain gas at z = 0, after UV ionisation, 
stellar feedback was the most prominet mechanism, while 
tidal and ram pressure stripping were close in magnitude. 
Finally, UV ionisation was the most important for the sub- 
halos that never formed stars, with some impact from ram 
pressure stripping. Tidal stripping and stellar feedback were 
negligible. Note the caveat: since it is impossible in all cases 
to clearly distinguish the mechanism that leads to the loss of 
a gas particle, the boundaries between the mechanisms are 
not clearly defined and the percentages should, therefore, be 
take as indicative of the relative importance of the various 
gas loss mechanisms. 

For the lower mass subhalos, the amount of mass lost 
adds up nearly to the cosmic baryon fraction (« 0.17), in 
part because our analysis relied on pairing dark and gas 
particles. However, not only do the higher mass subhalos 
contain more baryons than the cosmic fraction, they contain 
more stars than the cosmic fraction. To understand how the 
higher mass subhalos form so many stars, we investigated 
the origin of these stars. Figure [17] shows the mass evolution 
of a 7.1 X 10^ Mq subhalo that ends up with more than 
the cosmic baryon fraction in stars. The mass evolution is 
divided into categories based on whether the particles were 
twins of the dark matter present at the maximum mass. 
While most of the stars formed from gas that was a twin of 
this dark matter, almost 10% of the stars formed from gas 
that were twins of dark matter that were not members of 
this halo at its time of maximum mass, or any of the outputs 
immediately before and after the time of maximum mass. 

Figure [18] shows how this extra gas (marked as light 
green) comes from a much wider region than the dark mat- 
ter (marked as brown). While such accretion could be a nu- 
merical artifact of over-efficient gas cooling, it could also be 
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no gas but 


gas 






never stars 


have/had stars 
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Min Mass 
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1.10 X 10* 


7.79 X 10* 


3.14 X 10" 


Max Mass 


(Mq) 


2.09 X 10" 


4.30 X 10" 


6.03 X 10"' 


UV (%) 




15.73 


11.19 


6.33 


Ram (%) 




1.10 


1.54 


0.26 


Tides (%) 




0.12 


1.87 


1.93 


Sfb (%) 




0.01 


2.08 


0.64 



Table 1. The cumulative effect of each gas-loss mechanism (UV 
ionisation, ram pressure stripping, tidal stripping, and stellar 
feedback) given as an approximate percentage of the subhalos' 
maximum mass, averaged over all the subhalos in each category. 
The first category are the subhalos that had no gas at 2 = and 
never formed stars (represented by subhalo (c)), the second are 
those that have no gas at 2 = but had stars at some point 
in their history (represented by subhalo (b)), and the last cate- 
gory are those subhalos that retain both gas and stars at 2 = 
(represented by subhalo (a)). Also included are the minimum and 
maximum masses of the subhalos in each category. 
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Figure 17. Time evolution of a 7.1 X lO" Mq mass subhalo's mat- 
ter, broken down into dark matter present (twin) and not present 
(not twin) at the time of maximum mass, gas twinned or not 
from the maximum mass dark matter, and stars that did and did 
not come from twinned gas. The dashed vertical line is the time 
of maximum mass. Despite the difference in magnitudes between 
the twin and non-twin dark matter, the twin and non-twin gas 
particles are comparable, suggesting that the dark matter draws 
on gas outside its region of origin. 



a unique feature of satellites that orbit in high-density re- 
gions like a massive galaxy's hot halo. The mechanism that 
appears in the simulations is that high mass satellites quickly 
form stars out of gas that are twins of member dark matter 
particle. After they form stars and the stellar feedback cools 
down, there is less gas to provide pressure support to keep 
hot gas from the main halo out of the satellite. So this gas 
is accreted, cooled, and finally forms stars. 



5 CONCLUSIONS AND DISCUSSION 

To gain insight into the missing satellites problem, we com- 
pared the satellite luminosity functions of two simulated 




Figure 18. Snapshots at 2 = 6 of a 7.1 X lO" Mq mass subhalo 
(left) and a 2.4 X lO" Mq mass subhalo (right). The subhalos 
are in brown, while light green is the gas that will produce all 
the non-twin stars that will end up in the subhalos. This gas will 
eventually converge into the subhalos. The subhalo on the left 
will end up in a tight orbit around the host, and hence the several 
shells of gas that it will draw upon. The background colours are 
gas temperature as in Figure [T0I 



galaxies from the MUGS project (g5664 and gl5784) with 
late-type galaxies. 

The cumulative number of luminous satellites in gl5784 
was only slightly higher than that observed in the Milky 
Way, though there were an excess of high luminosity satel- 
lites that created a "knee" in the satellite luminosity func- 
tion that is not observed. Other SPH simulations of similar 
or lower resolution to ours have found that the missing satel- 
lites problem is no longer a matter of an order magnitude 
difference between the Local Group and simulated subhalo 
populations. When we compared our luminosity function of 
gl5784 to a simulation of the same galaxy at a lower gas res- 
olution both luminosity functions were similar in their area 
of overlap down to Mv ~ —8.2, though the low resolution 
run did suffer from having fewer gas particles than needed 
to properly resolve star formation and therefore it had fewer 
stars. A couple of our dwarfs had luminosities comparable 
to the recently discovered ultra-faint dwarfs, but at the res- 
olution of these simulations, they had only one or two star 
particles, so it is impossible to draw any conclusions about 
the formation of fainter dwarf galaxies from these simula- 
tions. 

The satellite mass function of gl5784 revealed a large 
population of dark satellites. In our more massive galaxy 
gl5784 (1.4 X 10^^ M0) the subhalos constituted 6.0% of 
the host galaxy's mass, while g5664's subhalos were 4.4% 
of the host galaxy ' s mas s. This fits within the range that 
iDalal fc Kochanekl (|2002l ) found from probing substructure 
with gravitational lensing, between 0.6% and 7.0%. 

We used two methods to determine how the dark satel- 
lites lost their baryons and the effect of negative feedback 
on the luminous satellites. 

One method was to simulate the less massive galaxy 
than the Milky Way, g5664, using several different physical 
treatments. The simplest included no UV or stellar feed- 
back. In the second, UV was added. These two simulations 
were compared with the standard MUGS simulation that in- 
cluded both UV and stellar feedback. The effect they had on 
the subhalo populations was significant. UV feedback alone 
stopped star formation in all the satellites with total masses 
less than 2 x l(f Mq. When stellar feedback was added, 
it reduced the luminosity of several additional subhalos so 



14 Nickerson et al. 



that only a couple of star particles formed in those subhalos 
before the feedback ejected all the remaining gas from the 
subhalos and eliminated the possibility of future star for- 
mation. In more massive subhalos, the stellar feedback had 
little impact on reducing the star formation efficiency. This 
unbalanced influence of the stellar feedback may have been 
due in part to the quantized feedback that was used in the 
MUGS simulations. 

The second method was to analyse the individual evo- 
lution of satellites in a more massive halo. We made a com- 
prehensive study of the mechanisms that remove matter 
from subhalos by defining criteria for mass loss due to the 
UV background, tidal stripping, ram pressure stripping, and 
stellar feedback. This analysis reiterated the strong impact 
ionisation had on low mass satellites. We used metals to 
track stellar feedback, and found that its impact was largest 
on subhalos of medium mass that had formed stars, with 
lesser impact on the highest mass subhalos, and no impact 
on the lower mass subhalos. 

A strange phenomenon was apparent in the higher mass 
subhalos. These subhalos contained a higher fraction of their 
maximum mass in stars than the cosmic baryon fraction. 
Subsequent analysis showed that accretion of baryons was 
not confined to the same limited region from which dark 
matter was accreted, but from a larger region surrounding 
the subhalo and even from across the hot gaseous host halo 
as the subhalo moved through its orbit. While this may par- 
tially be another symptom of overcooling that has been long 
noted in simulations, it may also point to the enhanced bary- 
onic accretion possible by subhalos in high density regions. 

Stripping, either ram pressure or tidal, also plays a vi- 
tal role in shaping the satellites that were analysed. Ram 
pressure removed whatever gas remained in small subha- 
los that had most of their gas removed during ionisation. 
In some cases, more gas was strippe d from subhalos than 
would have been predicted based on a lGunn k, GottI (| 19721 ) 
analysis. Tidal stripping removed most of the gas from the 
more massive satellites, making them comparable with Lo- 
cal Group dSphs rather than dirrs. The stripping became 
apparent once the subhalos crossed inside the virial radius. 

Tidal stripping was important for the subhalos that 
had the closest encounters with the main galaxy. In some 
cases, tidal stripping removed enough of the outer layers of 
dark matter that the total mass of the satellites dropped 
below the ionisation mass limit of 2 x 10^ ^0- Because 
tidal stripping reduces the total mass of subhalos by dif- 
ferent amounts, it is critical to organize the satellites by 
their maximum mass rather than their mass at 2 = to 
see a continuous behavior in the baryon fraction as a func- 
tion of mass. Many authors have not ed the similarity in 
mass inferred in Local Group dSphs |Bullock et al ] I2OOOI : 
IStrigari et al ] I2OO7I : IPeiiarrubia et al.l I2OO8I ) bv extrapolat- 
ing satellite total masses using NFW density profiles. Our 
simulations point out that because of tidal stripping, these 
satellites may no longer contain that much mass. However, 
those extrapolated masses may be similar to the maximum 
mass of the satellite, and the constant lower mass limit sug- 
gests a mass-dependent gas removal mechanism like ionisa- 
tion. What is not clear in the simulations is why the effi- 
ciency of star formation varies so much from dSph down to 
ultra faint galaxies if they did form from subhalos that were 
all the same mass. 



We should also extend this work to include all MUGS 
galaxies to determine if there are specific factors in the en- 
vironment or history of individual galaxies that affect the 
satellite population, as well as compare the subhalo popu- 
lations outside of the main halos to the ones that end up 
within the virial radius of their main halos. 
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